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The problem of deep laser cooling of ^“^Mg atoms is theoretically studied. We propose two-stage 
sub-Doppler cooling strategy using electro-dipole transition 3®P2—^-S^Ds (A = 383.9 nm). The first 
stage implies exploiting magneto-optical trap with and a~ light beams, while the second one uses 
a linLlin molasses. We focus on achieving large number of ultracold atoms (Te//<10 ^K) in a cold 
atomic cloud. The calculations have been done out of many widely used approximations and based 
on quantum treatment with taking full account of recoil effect. Steady-state average kinetic energies 
and linear momentum distributions of cold atoms are analyzed for various light field intensities and 
frequency detunings. The results of conducted quantum analysis have revealed noticeable differences 
from results of semiclassical approach based on the Fokker-Planck equation. At certain conditions 
the second cooling stage can provide sufficiently lower kinetic energies of atomic cloud as well as 
increased fraction of ultracold atoms than the first one. We hope that the obtained results can assist 
overcoming current experimental problems in deep cooling of ^"^Mg atoms by means of laser fields. 

Cold magnesium atoms, being cooled in large number down to several /rK, have certain interest, for 
example, in quantum metrology. 

PACS numbers: 37.10.De, 05.10.Gg, 06.30.Ft 


I. INTRODUCTION 

Laser cooling and trapping of neutral atoms plays 
important role for many directions of modern quantum 
physics. One of the directions is quantum metrology 
that experiences galloping progress nowadays. It is aimed 
on producing of standards of various physical quantities 
and on carrying out precise measurements with the help 
of them (e.g., see [l|). At present the most precise 
measurements are possible for the physical quantities 
as frequency and time. It is due to success achieved 
in producing etalons (standards) for these quantities. 
Contemporary time standard is based on a frequency 
standard, which defines its stability and accuracy to a 
considerable degree. At that, frequency etalons can be 
used not only as a basis for time standards, but also 
for conducting precise measurements of other physical 
quantities and constants as, for instance, electrical 
current and voltage, magnetic field, length, Rydberg and 
fine-structure constants. 

High-accuracy experiments for versatile examination 
of relativistic and quantum theories have become 
feasible owing to modern frequency standards. Among 
practical applications of time and frequency standards 
the broadband communication networks, navigational 
and global positioning systems should be mentioned 
especially. Many laboratories in word-known scientific 
centers do research in the field of frequency standards. 
One of the latest trends in this field is connected with the 
concept of intercity or even international quantum clock 
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network that could combine time etalons from various 
laboratories and countries into one system M- 

There are two main directions of primary frequency 
standards development: the technology of a single ion 
confined in a electro-quadrupole trap and the second 
one based on many neutral atoms trapped in an optical 
lattice (e.g., see [ 3 , 0). The latter direction is much 
newer than the former and it undergoes intense progress. 
The idea of neutral atoms trapping in a periodic light 
potential is not new and it was actively studied in 
1970s (see monograph Q and citations in it). As 
for metrological purposes this idea has experienced the 
second birth in the beginning of XXI century after 
noticeable progress in technique and methods of laser 
cooling of atoms, development of the “magic”-wavelength 
concept Eini, and also experimental and theoretical 
success in the field of spectroscopy of forbidden atomic 
transitions [l2l - [l^ . At present, stability of optical- 
lattice-based frequency standards is on the same level 
with single-ion standards and in some cases even better. 
The state-of-the-art prototypes reached instability and 
uncertainty on the relative levels of 10“^^—10' -18 

HI- 

Ones of the main candidates for producing the new- 
generation frequency standards are alkaline earth and 
alkaline-earth-like atoms: Yb (for instance, see 1211- - 
^), Ca m, Sr [H [H, Hg [H] and Mg [23, 
These atoms are the most appropriate because 
of narrow spectroscopic lines due to forbidden optical 
transitions from the ground state ^Sq to the lowest 
excited triplet state ^Po,i ,2 (see Fig. [T]). Moreover, one 
more key circumstance consists in the existence of so- 
called “magic” wavelength for these transitions. Under 
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the magic-wavelength optical field the linear (in the 
intensity) light shift is canceled. Also, one of the last 
tendency in this field is connected with spectroscopy of 
transition ^Sq^^Pq in even isotopes (with zero nucleus 
spin), which is highly forbidden. Frequency of this 
transition is immune to many frequency-shift effects, 
therefore this transition can be exploited as a good 
“clock” transition. In spite of the transition is highly 
forbidden, it has already been observed by means of 
magnetic-field-induced spectroscopy in [2lj, 

®®Sr [2i,[33-lll and ^^Mg H. 

To date, atoms of the first four elements (Yb, Ca, 
Sr and Hg) can be effectively cooled by the laser 
methods down to ultralow temperatures, approaching 
to the recoil energy limit, what is required for effective 
loading an optical lattice (~ 1-10 ^K). Besides, sub¬ 
recoil temperatures can be obtained by using evaporative 
coo ling technique, getting Bose-Einstein condensation 
[ 33 -I 33 . Unfortunately, researchers have not been able 
to reach the same great success with Mg atoms. In 
particular, neither two-photon laser cooling nor laser 
quenching [s^ methods have appeared to be ineffective 
in the case with magnesium. The minimum temperature 
of a magnesium cloud that has been obtained by laser 
cooling method is about 500 /iK, what is rather far 
from desirable range of values, in particular, from the 
recoil temperature (3-10 /rK, depending on an atomic 
transition). 

At the same time, magnesium atom has some 
advantages with respect to the other candidates for the 
frequency standard. Thus, black-body radiation (BBR) 
shift is one of the main limiting factors for accuracy 
and stability of quantum frequency standard (e.g., see 
BBR shift of the clock transition 
for magnesium is much smaller than for Yb, Ca, Sr 
and just a little bit higher than for a mercury atom 
(see Tab.l). However, from an experimental viewpoint 
Mg has some advantages against Hg. In particular, 
mercury atom requires noticeably smaller wavelengths 
for laser spectroscopy, cooling and trapping than Mg. 
For instance, effective trapping of cold atoms in a 
nondissipative optical lattice and producing Lamb-Dicke 
regime need for an optical potential depth at the level 
of 50-300 in the recoil energy units (e.g., see @,111). 
It means that a highly intensive laser field at the magic 
wavelength Am should be applied. But, as it can be seen 
from Table 1, Am (Mg) Ri 468 nm and Am(Hg)Ri363 nm. 
Therefore, from an experimental viewpoint, producing 
deep optical potential for mercury atom is more difficult 
problem than for magnesium one due to the much 
smaller wavelength. Besides of relatively small BBR 
shift, magnesium atom has one more advantage with 
respect to Ca and Sr. It consists in absence of optical 
pumping of atoms on the non-resonant level 3^D2 during 
the realization of laser precooling with the help of strong 
dipole transition 3^So—(see Fig. ^ to reach the 
temperature down to a few millikelvins |27l . [ 2 ^, 113 • 

Recent experiments [H, 113 showed some progress 
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FIG. 1. Partial energy diagram of ^^Mg atom. Solid 
lines denote the cooling transitions with corresponding 
temperatnre limits, while dashed lines denote possible “clock” 
transitions, which can be used for laser stabilizing. 


TABLE I. Data for several atomic elements relevant for new- 
generation freqnency standards: Xd is a wavelength of the 
clock transition 3^So—^-S^Po and Am is its magic wavelength, 
BBR frequency shifts are indicated with respect to absolute 
frequencies of clock transition {Xd is taken from NIST Atomic 
Spectra Databas^B)- 


Atom 

Ac/ 

Am 

BBR shift 

Sr 

698.5 

813.5 [IJ 

-5.5x10”^® 

Yb 

578.4 

759.4 [21,^ 

-2.6x10"^® 

Ca 

659.7 

735.5 [2^ 

-2.6x10"^® 

Mg 

457.7 

« 468 [40] 

- 3.9 X 10-^® [Tl] 

Hg 

265.6 

362.6 [2^ 

-2.4x10”^® [42. 43 ] 


^ http://physics.mst.gov/asd 


in cooling of ^^Mg atoms. The atoms were cooled 
down to the record temperature equaled to 1.3 /iK and 
confined in an optical lattice. However, the final number 
of atoms was of the order of 10"^, what was about 
0.01% from the initial number of atoms in a magneto¬ 
optical trap (MOT), involved the cyclic triplet dipole 
transition 3^P2^3^D3. That great loss in atomic number 
was due to the fact that velocity selection technique, 
having similarities with evaporative cooling, was used 
for reaching such ultralow temperature, but the laser 
cooling method in MOT, unfortunately, showed the cloud 
temperature equaled to only about 1 mK. That result 
noticeably yielded to successful results with the other 
atomic elements (Ca, Sr, Yb, Hg). At the same time, we 
believe that laser cooling strategy for magnesium atoms 
can be proper tuned for getting much better result of 
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laser-cooling temperature as well as number of atoms 
trapped. 

Therefore, we can state that problem of deep cooling 
of magnesium atoms by means of laser radiation is 
still unsolved. Moreover, increasing of ultracold atomic 
number has principal importance for many applications 
of cold atoms. For instance, authors of the paper 
managed to obtain Bose-Einstein condensation composed 
of strontium atoms. Besides, frequency-standard 

stability depends on an atomic number in an OTtical 
lattice and it increases with the number increases [I|, . 

All these things considered, we can conclude that it is 
important to solve the problem of deep laser cooling 
of magnesium atoms (down to T^l-10 fiK) as well as 
to provide much larger number of ultracold atoms in a 
lattice. 


II. LASER COOLING IN MOT: 

SEMICLASSICAL APPROXIMATION 

Laser cooling of neutral atoms in a magneto-optical 
trap is one of the main cooling methods. At first time the 
laser field composed of six beams with orthogonal circular 
polarizations {a^a~ configuration) was suggested in 
[49l| as effective way for simultaneous cooling and 
trapping of atoms. Narrow spectral lines allow laser 
cooling of various atoms down to a few tens and 
units of microkelvin, and even lower. In particular, 
narrow intercombination transition 4^So—in '^°Ca 
(2~27rx400Hz) provides temperatures around 4-6 /iK 
[50, im just by the help of Doppler cooling process. There 
were good results with intercombination transitions also 
for the other elements: Sr [s^, Yb and Hg [d^, for 

even as well as for odd isotopes. In certain aspects even 
isotopes, having a zero nuclear spin, are more attractive 
for frequency standards of new generation. However, as it 
has been already noted in the Introduction, still there are 
not satisfactory results of cooling ^'^Mg atoms by means 
of laser radiation in contrast to the other elements. 

The first our attempt to solve the problem with 
deep laser cooling of magnesium was undertook in the 
recent work [s^, where the detailed theoretical study 
of magnesium kinetics in ID MOT using the dipole 
transition 3 ^P 2 —^-S^Ds was conducted. The theory was 
based on the semiclassical approach iil , based on the 
well-known assumptions: 

Wr-ec < min{7,7S'} (1) 

and 

Ap 3> hk . (2) 

Here ijjrec=hk‘^/2M is the recoil frequency, M is mass 
of an atom, k=2'!TlX is wave number. The saturation 
parameter S is defined as 


^ (7/2)2 + ^2 ’ (3) 

where 7 is the spontaneous relaxation rate of excited 
state, 6=uj—uJo is the detuning of laser radiation 
frequency w from the transition frequency ojq, and H is 
the Rabi frequency. 

Condition o implies that recoil frequency must 
be rather small in comparison with a typical rate of 
establishment of steady state among atomic internal 
degrees of freedom. In particular, in the case of an atom 
without any degeneracy of the ground state this rate is 
defined by 7 . If there is a degenerate ground state and 
optical pumping can occurs, this rate is defined by 7 or 
the pumping rate jS, depending on what is smaller. The 
second semiclassical requirement ([ 2 ]) implies that typical 
width of stationary linear momentum distribution f(jj) 
must be much larger than the recoil momentum from 
emmision/absorption of a photon. 

Doppler limit for temperature of laser cooling T/j, 
which can be achieved at the frequency detuning 
i5=— 7 / 2 , can be figured out from equation for minimum 
kinetic energy in one-dimensional case: 


Tpmin 7 rri 

^kin — — 


40 


hr). 


(4) 


Strictly speaking, this equation is valid for transition 
Jg=0^Je=l. It was found in under a'^a~ 

configuration (also see da). If we use this formula 
for getting estimate of Td in the case of transition 
3 ^P 2 —i> 3 ^D 3 (7ft!27r26.7MHz), we immediately find 
425 /iK. For effective trapping of atoms with such 
relatively high temperature the large intensity of cw 
optical lattice field at the level of tens of MW/cm^ 
is required, what is hardly feasible in an experiment. 
Therefore, much lower temperature of atomic cloud 
is needed. At the same time, since the transition 
considered has degenerate energy levels, one can 
anticipate activation of so-called sub-Doppler mechanism 
during the laser cooling in MOT under the polarization- 
gradient field. In principal, this process would overcome 
the Doppler limit (|3]) and show much lower temperature 
than in the case of Jg= 0 —>'Je=l. 

Semiclassical approach is based on kinetic equation of 
Fokker-Planck type on the Wigner distribution function 
in phase space f{z,p). That equation can be acquired 
by reducing of exact quantum kinetic equation on 
the density matrix in the series on small parameter 
hk/Ap<^l until second-order terms. This procedure is 
well-known and it has been done by many authors (e.g., 
see Eventually the following equation can be 

obtained: 


[-^Fiz,p) + ^Diz,p)\ f{z,p). (5) 
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FIG. 2. Momentum distributions of magnesium atoms 
at (5= —57~—27rx 130 MHz, 7=20mW/cm^ (solid) and 
7=470mW/cm^ (dashed). 


Here F{z,p) is laser-field force on an atom, D{z,p) is 
diffusion of an atom in the light field. This equation 
must be completed with normalizing condition that in 
one-dimensional periodic laser field has the form; 


1 

A 


-t-A/2 +00 


dz 


f{p, z)dp = 1 . 


-A/2 -oo 


One-dimensional cr+cr" laser-field configuration allows 
significant simplifying of ([5]), at that the dependence / 
on z vanishes (see section ITlI Bl) . 

Our semiclassical calculations have been done 
beyond many widely used approximations (for instance, 
slow atoms and weak field approximations). As it 
has been shown the minimum kinetic energy achievable 
in MOT is close to SOxErec, where Erec=^rec is 
the recoil energy. The effective temperature, which 
can be associated to this value, TeZ/wlSO/rK. It is 
approximately three times lower than the estimate of 
Doppler limit 425 /iK, but, unfortunately, it is still 
very far from desirable range of values and, in particular, 
the recoil temperature Trec=5/rK. 

Let us consider a question on validity of semiclassical 
approach to the magnesium problem. Indeed, as it will 
be shown further on the basis of semiclassical treatment, 
the optimal parameters of cooling field for transition 
3 ^P 2 ^ 3 ^D 3 can be chosen as i5=—27rxl30 MHz 
and 7=500 mW/cm^. The corresponding saturation 
parameter S'r;4xI 0“^. In spite of such low saturation, 
the first semiclassical requirement m is still satisfied, 
because a;rec= 2 xl 0 “^ 7 . At the same time, typical 
momentum distribution width Ap may not satisfied the 
second semiclassical condition ©• Indeed, in general 
case the distribution can have complex shape. In 
particular. Fig. [5] shows two examples of momentum 
distributions for different values of light field intensity. 
The distribution acquires two-peaked profile at low 


intensity 7=20mW/cm^: there is the high-contrast spike 
on top of the wide background. This background 
conditionally describes “hot” fraction of atoms in a cloud 
with effective temperature Tgf1—10 mK, while the 
spike corresponds to ultracold fraction with T^ff'^ 1 /rK. 
Similar distributions were observed earlier (e.g., see 621 
with semiclassical low-saturation-limited calculations of 
Sisyphus cooling of Cs atoms with the help of transition 
F"g=4—>-F).=5 or the quantum-treatment calculations for 
atomic W-type scheme in (^). In our case the narrow- 
spike width is about hk and the requirement Apl^hk 
is not satisfied at all. With increasing the intensity 
(7=470mW/cm^) the two-peaked shape disappears. 
However, the distribution as a whole is still sufficiently 
narrow and the second requirement @ is satisfied with a 
good margin. Also it should be noted that the condition 
m, as a matter of fact, depends on the value of total 
angular momentum Eg. In other words, at the same 
saturation parameter S the requirement ([2]) can be valid 
for small Eg and getting not valid with its increasing. 

Basing on the aforesaid, we can conclude that more 
precise theoretical treatment is needed in the case 
with magnesium for adequate description of kinetics of 
ultracold atoms. This treatment can be based on the 
density matrix formalism with full account for the recoil 


effect (e.g., see 


M)- Moreover, as it will be 
seen in the next section, the quantum-treatment results 
noticeably differs from the semiclassical ones, based on 
the equation ([5]). That difference, in particular, gave us 
an idea for exploiting the second stage of sub-Doppler 
laser cooling for getting the desirable results. 


III. FULL ACCOUNT FOR THE RECOIL 
EFFECT 

Let us consider the problem of laser cooling of 
magnesium atoms out of semiclassical approximation 
limit as well as some other widely used approximations 
(weak-saturation limit, secular approximation, etc.). 


A. Problem statement 

We assume the laser field to be one¬ 
dimensional, composed of two plane monochromatic 
counterpropagating light waves with equal frequencies 
and amplitudes (the quantization axis z is collinear to 
the wave vectors): 


E(z,t)= -H F;oe 2 e-*(‘^*+'=^) -h c.c. = 

= Eoe{z)e~^‘^^ + C.C., ( 6 ) 

where ei _2 are the unit complex vectors of waves’ 
polarizations, while e(z) is the following complex vector 


e(z) = eie*'=^ -f 626 "*'=^ . 


(7) 
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Nonzero components of the vectors ei ^2 in the spherical 
basis are 


= — sin(ei — 7r/4), = — cos(£i — 7r/4), 

62^ = — sin(£:2 — 7r/4) 

= — cos(e 2 — 7r/4) (8) 

Here £12 are the ellipticity parameters (in particular, 
£=±7r/4 corresponds ro righ- or left-circular polarized 
wave, e=0 is for linear polarization), <p is the angle 
between main axes of polarization ellipses. For instance, 
the case with £ 1 . 2=0 and (^=7r/2 corresponds to linJ^lin 
field configuration. 

Here quantum treatment of atomic kinetics under the 
laser field (l6|) is based on the equation on single-atom 
density matrix in coordinate two-point representation 
that has the form (e.g., see [1, [H, Q): 


dp{zi,Z2,t) 

dt 


hi 


H{zi,t)p-pH{z2,t) +f{p}, (9) 


with the Hamiltonian 


H{zi,t) = {pf/2M) + Ho + V{z,,t). (10) 

The first term in the Hamiltonian is the operator 
of kinetic energy of an atom {pi is the linear 
momentum operator), Hq describes intratomic degrees of 
freedom, operator V corresponds to the atom-field dipole 
interaction, and the linear operator functional r{. .. } is 
respective for relaxation processes in an atom. Let us 
introduce the projection operator onto the excited atom 
state: 

P"" = ^^1 Fe,rne){Fe,me \ , ( 11 ) 

me 


and the Wigner vector operator T, whose spherical 
components are: 




Fe , {Fg , THg 


( 12 ) 


me,mg 


with cr=0, ±1 and Cp'””)'’,, the Clebsch-Gordan 

■C g , ''I'g , J-O 

coefficients (e.g., see [H). Then the terms iJo and V 
from (US in the resonant approximation can be written 
as 


% = -h6P^ 


(13) 


and 


V{zi) = -hRT-e{zi) +h.c. = -hRV^o(^Zi) +h.c. (14) 


Here R=EQd/h is the Rabi frequency (with d the 
reduced matrix element of dipole operator of an atom), 
V’^^{zi)=T-e{zi) is the dimensionless operator of atom- 
field interaction, depending on the coordinate in general 
case, h.c. means Hermitian-conjugate term. 

Introduce the new coordinates: 


z = k{zi P Z2)/2, q = k{zi - Z 2 ) , (15) 

in which the spontaneous relaxation operator from ([2]) 
acquires the form: 

r = -|(p'=P + pP'=)+7 E Ud)TlpTa, (16) 

(T= 0,±1 


with 


2 V q 



sin(( 7 ) 

cos(g) \ 

) • 


COs{q)\ 

’ 


(17) 


Note that in the absence of recoil effect, i.e. in the limit 
q ^ 0, we have = 1. 

The density matrix can be divided into four matrix 
blocks: 


P = 


pgg \ 

^g j ■ 


(18) 


Matrix blocks and describes populations of the 
ground and the excited states as well as low-frequency 
(Zeeman) coherences. Blocks and pP^ are responsible 
for optical coherences. For the new coordinates m and 
using all introduced notations the new equations on the 
density matrix blocks can be easily acquired from ([^ . 
So, in the steady state we have 


-2iuJr-^-^^^{z,q) ^ Ca{q)f^pTa + 

+ii?[F"sl(z + |);3"s - ^"F"s(z - |)] , (19) 

N. 




= iR 


dqdz 

i/"9(z-K |)^«-^SF"9t(z- |)] , ( 20 ) 

52 V 




= iR 


i/"9l(z + |)p""-p®9F"sl(z-|)J , (21) 




= iR 


dqdzy 

+ |)p®9 - - |)j 


( 22 ) 
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These equations compose a basis for further theoretical 
analysis. For instance, probability density of atoms in 
the momentum space can be found from the formula: 


+ 00 TT 

fip) = j dqj dzTr{p(z,g)}e-*J>« . (23) 

— OO —TT 

Here the linear momentum of an atom evaluated in 
the recoil momentum units hk and Tr[...] denotes 
trace operation. Momentum distribution f{p) must be 
normalized: 


FIG. 3. Transformation the old coordinate frame K to the 
new one K'. 



+ 00 

j f{p)dp = 1. (24) 

— OC 

This means that the set of equations (IT^ - (l2^ must be 
supplemented with the condition: 

TT 

^ / Tr{p(z,Q = 0)}dz = 1. (25) 

— TT 

Average kinetic energy of an atom in the recoil energy 
units can be evaluated, for instance, with the help of the 
following formula: 


The considered field configuration has some unique 
features. First of all, the field has homogeneous intensity 
(it does not depend on z-coordinate). And the second, 
the field polarization also can be made homogeneous 
(e.g., see [H, H^). Indeed, let us pass to the new 
coordinate system K', in which z'-axis coincides with 
z-axis in the old AT-system, while the axes x' and y' 
rotate around z-axis by the angle q;=— z=—A:(zi-|-Z2)/2 
(see Fig. [3]). In the AT'-system the linearly polarized 
total-field vector does not rotate anymore (without loss 
of generality it can be considered to be directed along x'- 
axis). Then in the new system the interaction operator 
F®® from (fT31) does not depend on the coordinate zp. 


+ 00 

Afc = / p^f{p)dp. (26) 


B. The first stage: Cooling in MOT 

As a rule, in magneto-optical trap atoms are localized 
in weak-magnetic-field region (in the vicinity of trap’s 
center). Therefore, magnetic field does not affect 
significantly on the temperature of a cloud and we 
omit it here from our analysis. In other words, we 
consider kinetics of atoms in ID laser field, composed of 
two counterpropagating beams with orthogonal circular 
polarizations (cr^cr” configuration). Then in equation 
dSl) we can take £i=7r/4, £ 2 =—7r/4, ip=Q and polarization 
vector of total light field in spherical basis takes the form: 


F"9(zi) ^ Vx[q) = , (29) 

^%{q) = n{q)- (30) 

Since the relaxation operator F from m also does not 
depend on z, the density matrix in the new coordinate 
system is not the function of z and it depends only on 
g-coordinate. This circumstance significantly simplifies 
numerical evaluations of the density matrix equations. 

The operator of rotation D{n, a) can be exploited for 
getting the equations on density matrix in the new basis 
(e.g., see [i^). Here the unit vector n defines a rotation 
axis, while a is a rotation angle. In our case it is natural 
to coincide n with quantization axis z and take a=—z. 
Then action of the rotation operator on a wave function 
I Fa,ma,Zi) reduces to a simple multiplication by 6*"*“, 
i.e. 


e(z,) = . (27) 

This form corresponds to the laser field with linear 
polarization, which rotates by an angle a=—kzi during 
propagation along z-axis. At that the dimensionless 
operator F®® from da is 


D{n, a) I Fa, nia, Zi) = 6*"*““ | Fa, m-a, Zi) , (31) 

with {a=e,g). In the system K' the set of equations 
(HU)-© take the form: 


2uJr 


dq - 


Fz,p{<l) =r{p(g)} + i5 




Vi{q)p- pV 2 {q) 


(28) 


-I- iR 


(32) 
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FIG. 4. Comparison the results of semiclassical (dashed) and qnantnm (solid) treatments at Sqss—27rx 130MHz. (a) 
Average kinetic energy of an atom as the function of light field intensity, (b) Ultracold fraction of atoms in a cloud. 



FIG. 5. Momentum distributions of magnesium atoms: 
comparison of semiclassical (dashed) and quantum (solid) 
treatments, (5=— 67 , i?~ 0.227 (/«20mW/cm^). 

At that the normalizing condition (l25l) becomes rather 
simple: 

Tr{p (<7 = 0)} = 1. (33) 

Figure |4ji shows average kinetic energy of an atom 
as the function of light field intensity, calculated 
on the basis of numerical solving the equation (1321) . 
Analogical dependence is also presented, gained by 
applying semiclassical approach on the basis of Fokker- 
Planck equation ([5]). As it is seen from the figure the 
semiclassical approach (dashed line) gives the minimum 
kinetic energy of an atom at the level of Emin^ ^OErec, 
what is several times smaller than the Doppler limit 
EuKi 87.5Erec- At the same time the quantum approach 
(solid line) shows the result for energy just a little 


bit smaller than the Doppler limit {Emin^G^Erec)- 
Hence the quantum treatment of the problem shows 
that it is hardly possible to cool magnesium atoms 
in MOT down to desirable range of temperatures on 
the basis of transition 3^P2—^■3^D3. All this agrees 
with the experiments of research group from the 
University of Hannover |4^. Effective temperature, 
corresponding to the minimum at the plot E{I) for 
quantum-treatment result, is about 310 fiK at frequency 
detuning d=—57«—27rx 130 MHz and light field intensity 
iRsllOOmW/cm^. 


Beside the temperature of an atomic ensemble it is also 
important to know a profile of momentum distribution 
of atoms in a cloud. It may be found very useful, in 
particular, for realization of evaporative cooling stage 
for achieving ultralow temperatures (~ 1 fiK). Let us 
consider a group of atom in the momentum space with 
p < 3 hk. Tentatively speaking we call this fraction as 
“ultracold” one. Figure shows number of atoms in 
the ultracold fraction Nc as the function of light field 
intensity I. As it is seen from the figure there is a 
maximum in vicinity of 500 mW/cm^. It should be 
noted that position of this optimum is not immediately 
the same as for the minimum of the dependence Ekinil). 
FigureHJo demonstrates that about 40 % of atoms can be 
concentrated in the ultracold fraction. For comparison 
analogical dependence is presented, calculated on the 
basis of semiclassical approach (dashed line), which 
lies noticeably higher than the former one. The 
dependencies Ekin{I) and Nc{I) lose to the semiclassical 
ones, because quantum treatment provides significantly 
different result for the momentum distribution in the 
vicinity of p«0. Indeed, Fig. [S] shows very sharp spike 
in the semiclassical case and a tiny peak as the result of 
quantum calculations. 











FIG. 6 . (a) Average kinetic energy of an atom under linA-lin light field configuration calculated on the basis of quantum 
treatment, (b) ultracold fraction of atoms in a cloud. Light field detunings are 5= —27 (dashed) and 5=—5'y (solid). 


C. The second stage: Cooling in an optical 
molasses 

Fortunately, solution of the problem of deep laser 
cooling of magnesium atoms can be found by involving 
the second stage of sub-Doppler cooling with the help 
of one-dimensional optical molasses. The molasses is 
composed of two counterpropagating light waves with 
orthogonal linear polarizations (lm_Llm configuration). 

In contrast to a~^a~ configuration, in the case of 
linJ-lin field the total-field polarization transforms from 
linear to circular (and back) along z-axis (e.g., see (b^). 
Consequently, there is no any rotating transformation 
of coordinate frame K that would make density matrix 
independent of z-coordinate. Therefore, we must solve 
the set of equations (mu-dill) on matrix p{z,q). We 
have solved the equation numerically on the basis of 
matrix continued fractions method. The details of the 
method can be found, for example, in and we do not 
reproduce it here. Instead of that, we just present the 
numerical results. 

Figure [6^ demonstrates much lower minimum kinetic 
energy than in the case of cr+cr“ field (see Fig. 0^, 
solid line). In particular, the minimum corresponds to 
Epil&Erec at JssSOO mW/cm^ (Te//~ 80/rK). Besides, 
as it is seen from Fig. [6 }d, ultracold fraction of atoms 
under linElin light field can be higher than in the case 
of a'^a~ (compare with Fig. IDo). Narrow structure 
in momentum profile in vicinity of p~ 0 becomes more 
visible than under cr+cr” field (compare Fig. [7] and 
Fig. [5]). Therefore, the second cooling stage involving 
optical molasses can provide lower temperature as well 
as larger number of atoms in ultracold fraction (up 
to 60%). After the second sub-Doppler cooling stage 
atoms may be loaded, for instance, to a dipole trap. 
At that the “hot” fraction of atoms (wide background 
at Fig. [7]) can be evaporated by proper choice of the 



FIG. 7. Quantum calculations of momentum distributions 
at (i=— 57 . Field strengths: 7 ?« 0 . 227 , 7«20mW/cm^ (solid) 
and i?~1.137, 7«600mW/cm^ (dashed). 

light potential depth, saving only the ultracold fraction 
in a trap (with effective temperature ~ 1 /iK). It should 
be noted that realization of the second sub-Doppler 
stage should eventually provide much more number of 
ultracold atoms in a dipole trap (or an optical lattice) 
after evaporation than without this stage. 


IV. CONCLUSION 

In conclusion we would like to summarize main results 
of the work. We have suggested using the second sub- 
Doppler cooling stage to solve the problem of deep 
laser cooling of magnesium atoms. The first stage 
implies using of magneto-optical trap involving dipole 
transition between triplet states 3^P2 and S^Ds. In 
particular, this stage was used in the experiment of 
researches from Hannover University In spite 
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of that the level 3^P2 is degenerate and one could 
anticipate activation of effective sub-Doppler mechanism 
of cooling in polarization-gradient field however, 

the conducted theoretical analysis has figured out the 
minimum temperature at the level of 310 //K, what is 
just a little bit lower than the estimate for Doppler limit 
of cooling (TDRi425^K). To reduce this value by several 
times we have proposed using the second laser-cooling 
stage with the help of optical molasses composed of 
two counterpropagating orthogonally linearly polarized 
waves (linJ-lin field configuration). In contrast to a^a~ 
field, applied in MOT, the optical molasses demonstrates 
much lower temperature (80 fJ-K). At the same time, 
the minimum achievable temperature of laser cooling 
in the case of ^“^Mg still is noticeably higher than 
for some other atoms, where sub-Doppler mechanism 
provides much better results. Most likely, it is due to 
relatively large recoil energy of magnesium atom. For 
instance, for considered transition in magnesium recoil 
frequency Wrec~ 27rx53kHz, while for ^^^Cs (^Si/ 2 , F=4 
—>■ ^P 3 / 2 , F=5) this frequency is significantly smaller 
(cur-p.r^27rx2kHz), what allowed cooling cesium atoms 
down to 2.5 

Besides temperature (average kinetic energy) of 
ensemble of atoms we have also paid attention to 
the linear momentum distributions in the steady state. 
In particular, we have investigated the problem of 
increasing concentration of atoms in ultracold fraction 
(a region in momentum space in the vicinity of p=0). 
Conducted numerical calculations have revealed the 
optimum parameters of laser field for maximization of 
the ultracold fraction (Te//~1 /iK). This fraction can be 
easily localized in a optical trap, while the other fraction 
(“hot” atoms) can be evaporated from the trap by proper 
choice of the optical depth. At that, it is the second stage 
of sub-Doppler cooling that can provide great increase 
of ultracold atomic number in comparison with the case 
with only the first stage realized (as in the experiments 

mm)- 

In our theoretical analysis quantum treatment with 
full account for the recoil effect has been exploited, i.e. 
we have not been limited by semiclassical or secular 
approximations as well as weak-field limit. It has 
allowed us studying kinetics of cold magnesium in a wide 
range of intensity and frequency detuning to determine 
the optimum parameters of laser field. Also we have 
compared data provided by quantum and semiclassical 
approaches. As the result of that comparison we can 
conclude that semiclassical approach in the case of 


transition 3^P2—^-d^Da in ^^Mg is not valid for adequate 
understanding the kinetics of ultracold magnesium atoms 
for a wide range of light-held parameters. Moreover, we 
can also conclude that for getting the adequate estimate 
of cooling parameters and understanding the problems in 
deep laser cooling of atoms it is quite necessary to treat 
the problems with the help of quantum approach. 

At the end we should note that in spite of theoretical 
analysis has been done out of many widely used 
approximations, we have assumed the problem to be 
one-dimensional. However, light-held conhguration used 
in a magneto-optical trap is always three-dimensional 
(three pairs of circularly polarized beams). Therefore, 
obviously, the results of such ID analysis may differ 
from the real experiment with 3D held. For example, 
one can refer to the papers [H, for getting the 
estimate of such kind of difference. In these papers 
calculations were done for ID and 3D conhgurations 
by the example of simple transition Fg=0^Fe=l in 
limits of semiclassical and slow-atoms approximations. 
At the same time, an optical molasses (the second 
cooling stage suggested), which is of the most interest 
from the viewpoint of deep laser cooling, can be 
implemented in 3D as well as in ID conhguration. Three- 
dimensional optical molasses for various transition of 
the type Fg=F^Ff.=F+l was investigated in ref. [t^ 
under weak-saturation approximation and with the help 
of adiabatic reduction of density matrix equations to 
ground state. Unfortunately that approximation gives 
good results not for wide range of parameters S and 
I that can have an interest from the laser cooling 
view of point (e.g., see the work [^ . where results of 
adiabatic approximation were compared with results of 
full quantum treatment). Three-dimensional quantum 
treatment with full account for the recoil effect and 
beyond the aforesaid approximations is quite difficult 
task that requires separate study. 
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